Interplay between SDW and induced local moments in URu 2 Si2 
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Theoretical model for magnetic ordering in the heavy-fermion metal URu2Si2 is suggested. The 
17.5 K transition in this material is ascribed to formation of a spin-density wave, which develops due 
to a partial nesting between electron and hole parts of the Fermi surface and has a negligibly small 
form-factor. Staggered field in the SDW state induces tiny antiferromagnetic order in the subsystem 
of localized singlet-singlet levels. Unlike the other models our scenario is based on coexistence of 
two orderings with the same antiferromagnetic dipole symmetry. The topology of the pressure phase 
diagram for such a two order parameter model is studied in the framework of the Landau theory. 
The field dependences of the staggered magnetization and the magnon gap are derived from the 
microscopic theory and found to be in good quantitative agreement with experiment. 
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I. INTRODUCTION 

URu2Si2 is one of the most intriguing heavy-fermion 
compounds. It exhibits a sharp second-order transition 
at T m — 17.5 K, which has pronounced effect on thermo- 
dynamic and kinetic properties^'^^ though the neu- 
tron diffraction experimentalise and the X-ray mag- 
netic scattering measurements^ have produced evidence 
for only tiny staggered moments fj, ~ 0.02-0.04/is at 
Q = (1,0,0). Magnetic excitations observed by inelastic 
neutron scattering experiments^ are reasonably well ex- 
plained within the model with exchange interaction in a 
singlet-singlet Van-Vleck paramagnet. 11 This model fails, 
however, to give a consistent description of the ordering 
temperature and small ordered moments unless the ex- 
change interaction J7q is accidentally close to a critical 
value (see section IV below). A weak antiferromagnetic 
ordering of a Van Vleck paramagnet cannot produce a 
measured jump in the specific heal£ and an electrical 
resistivity anomaly^ These experimental features rather 
resemble formation of a spin-density wave (SDW) , which 
involves approximately a half of the Fermi surface. In its 
turn, the SDW scenario is inconsistent with a longitudi- 
nal polarization of the magnetic excitations. As a result 
of this contradiction, various theoretical scenarios have 
been proposed in the past to explain intriguing behavior 
of '.'Pi.- t 1113,14 - 15,16,17,18,19,20,2122,23 

The proposed theories of URu 2 Si2 can be crudely di- 
vided into three broad categories. Phenomenological 
models of the first groupi^2£ take pragmatic approach 
and introduce a yet unknown primary or hidden order 
parameter ip, which drives the 17.5 K transition and is re- 
sponsible for a large specific heat anomaly at T m . Small 
ordered moments observed in the neutron experiments 
are described by a secondary Ising-like antiferromagnetic 
order parameter m. Depending on the symmetry of the 
hidden order parameter different coupling terms of the 
type ip m an d il^rr? are possible in the Landau energy 
functional. In particular, the model with a bilinear term 
predicts a nontrivial field-dependence of weak antifer- 
romagnetic moments with an inflection points Subse- 



quent neutron measurements have nicely confirmed such 
a prediction^ 

Second group of theories consists of specific proposals 
for the hidden order parameter. These models are fur- 
ther subdivided into two subgroups. The first subgroup 
includes models, in which the primary order parame- 
ter breaks the time-reversal symmetry. These are spin- 
density waves in higher angular momentum channels ^iiS 
triple spin correlatorsjiS^ orbital antiferromagnetism, 21 
and ordering of octupolar moments on uranium sites^ 
For such theories the bilinear coupling term with the an- 
tiferromagnetic vector is, in principle, possible, though 
some other crystal symmetries are required to be absent. 
The second subgroup includes theoretical models with a 
time-reversal invariant hidden order parameter such as 
quadrupolar— or spin nematic ordering^ In such a case 
the bilinear term is prohibited by symmetry and only a 
biquadratic interaction exists between the two order pa- 
rameters. 

The last third group of theoretical works includes at- 
tempts to realistically describe the microscopic inter- 
actions in URu2Si2. A so called Ising-Kondo lattice 
model 1 ^ describes interaction between conduction elec- 
trons in nested bands and local crystal-field split mo- 
ments. The mean-field calculation produces both a weak 
moment and an appropriate value of the transition tem- 
perature but does not reproduce the large specific heat 
jump. A closely related dual model— considers a sub- 
system of localized singlet-singlet levels and a subsystem 
of itinerant electrons with a similar assumption on the 
nesting condition. A better description of the thermody- 
namic properties of URu 2 Si2 has been achieved by adding 
the superexchange and the RKKY interactions between 
local moments. The field behavior still remains largely 
inconsistent with the experimental measurements^ Fur- 
ther development of the dual model has been recently 
suggested in Ref. |23- Note, that the above two micro- 
scopic models (i) completely neglect the Coulomb inter- 
action between charge carries and (ii) operate with a sin- 
gle (antiferromagnetic) order parameter. 

Up to now experiments have been unable to distin- 
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guish between the competing theoretical proposals. Two 
potentially important experimental developments pub- 
lished recently include pressure experiments and the 
NMR measurements. Investigations under hydrostatic 
pressure2£i2§i2! have shown that the P-T phase diagram 
is divided into two regions: a small moment antiferro- 
magnetic phase (SMAF) at low pressures and a large 
moment antiferromagnetic phase (LMAF) at high pres- 
sures with a first-order transition line in between. Such a 
discovery is generally consistent with phenomenological 
two order parameter scenarios for URu2Si2 (see section 
III, below) . Another experimental finding came out from 
the NMR measurements. Matsuda and co-workersS&Ssi 
have found that a paramagnetic Si 29 NMR absorption 
line persists well below T m accompanied by two much 
smaller peaks symmetrically shifted by antiferromagnetic 
field. Such an observation points at an inhomogeneous 
para-antiferromagnetic state below the 17.5 K transition. 
The peak intensities suggest that about 97% of the sam- 
ple volume is in a paramagnetic state. The puzzle of 
small uranium moments in URu2Si2 seems to be rein- 
terpreted as due to a phase separation between nonmag- 
netic state with a hidden order and small antiferromag- 
netic droplets with ordinary (large) value of staggered 
magnetization. Such a simple explanation of the main 
mystery of URu2Si2 is quite appealing but leaves with- 
out answer question why the inhomogeneous phase ex- 
ists not only in the vicinity of the first-order transition 
line at high pressures, but in the whole region of SMAF. 
The development of the antiferromagnetic Bragg peaks 
right below T m = 17.5 K seems to be highly acciden- 
tal in the phase separation scenario. Also, an inflection 
point in the field dependence of antiferromagnetic Bragg 
peaks cannot be explained if small peaks are purely due 
to a volume effect. Thus, physical implications of the 
NMR measurements^ 2 *! are not completely straightfor- 
ward and have to be further clarified. 

In this paper, we present a semimicroscopic model for 
URu2Si2, which is closely related to the above mentioned 
dual modelsfiii^ but operates with two order parameters 
in the spirit of phenomenological scenarios We also 
consider two magnetic subsystems: (i) local crystal-field 
split moments on U 4+ sites and (ii) conduction electrons 
in nested bands. In contrast to the previous worksiLi^ 
we assume that the electron-electron interaction is non- 
negligible and that it drives a SDW transition in the 
nested parts of the Fermi surface. The critical tempera- 
ture T m = 17.5 K is associated with Tsdw an d the SDW 
amplitude ip plays the role of a primary (hidden) order 
parameter. According to the LSDA calculations^ the 
nesting wave- vector is commensurate and corresponds to 
the experimentally observed two-sublattice antiferromag- 
netic structure. The SDW formed in conduction bands is 
responsible for large changes in thermodynamic and ki- 
netic properties of URu2Si2- At the same time, we argue 
that the SDW has a small form-factor and does not create 
significant Bragg reflection. Local polarization of ura- 
nium sites by a staggered magnetic field from the SDW 



induces tiny antiferromagnetic moments. The magnetic 
dynamics probed by neutrons is also determined by a 
localized subsystem. 

By making a single assumption about a hidden order 
parameter with the same symmetry as the observed an- 
tiferromagnetic ordering, we have further derived several 
results, which allow detailed comparison with available 
measurements and suggest future experimental tests: 

(i) The P-T phase diagram with line of the first order 
type transition, which terminates at the critical point 
below T m (P); 

(ii) Field dependence of staggered magnetization coin- 
cides perfectly with the experimental data for a realistic 
set of microscopic parameters; 

(iii) Field dependence of the excitation spectrum. 
The article is organized as follows. In section II, we 

introduce the model and discuss various features of the 
spin-density wave instability specific for URu2Si2. In sec- 
tion III, using the phenomenological approach with an 
appropriate Landau energy functional for two coexistent 
order parameters we investigate the topology of the P-T 
phase diagram of URu2Si2- In the next section, a system 
of localized crystal-field split singlet-singlet levels is con- 
sidered under combined influence of a uniform external 
magnetic field and a staggered "internal" field induced by 
the SDW. We calculate the field dependences of the stag- 
gered magnetization and the excitation spectrum. Com- 
parison to the experimental data is presented in section 
V, which is followed by discussion and conclusions. 



II. SPIN-DENSITY WAVE 

An unusual assumption made by phenomenological 
theorie a 15 ! 20 is a description of URu2Si2 with two or- 
der parameters having the same symmetry. On the first 
site such an assumption contradicts to a general spirit of 
the Landau theory of phase transition. Recent investiga- 
tions of a two-gap superconductor MgB2 have, however, 
demonstrated usefulness of the description of the super- 
conducting state by means of two weakly interacting s- 
wave condensates of the Cooper pairs (see, for instance, 
Ref . l3lh . The necessary condition for this is a significant 
mismatch of the pairing interactions in the two bands. In 
the absence of intcrband scattering of the Cooper pairs, 
each band has its own superconducting transition tem- 
perature. An interband interaction is always present in 
real metals and leads to a single transition into a state 
with two different gaps. The two gaps (order parameters) 
still keep different dependences on temperature, pressure 
and/or applied magnetic field. In relation to URu2Si2, 
the two order parameters ip and m should correspond to 
two significantly different magnetic subsystems. We sug- 
gest here that the primary order parameter may be an 
ordinary SDW. The common objections against a SDW 
transition are (i) smallness of ordered moments and (ii) 
longitudinal polarization of sharp magnetic excitations. 
This two properties can be reconciled with a SDW see- 
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nario by taking into account specific features of URu2Si2- 
In this section we discuss the former feature, while the 
properties of magnetic excitations are considered in sec- 
tion IV. 

The early experimental works on the specific heali 3 * 3 ^ 
and the magnetoresistance* in URu2Si2 have found 
strong evidences in favor of a charge or a spin-density 
wave instabilities in the heavy-electron subsystem at 
T m = 17.5 K. The fit of the electronic specific heat below 
the transition indicates that a gap A ~ 130 K devel- 
ops on 40% of the Fermi surface at T — * 0. This con- 
clusion has received a strong support from the de Haas 
van Alphen (dHvA) measurements^ Comparison of the 
measured dHvA frequencies to the ab-initio band struc- 
ture shows that two large pieces of the Fermi surfaces, 
band-18 hole and band-19 electron, are not observed 
at low temperatures, probably due to their partial re- 
moval below the ordering temperature. The above two 
sheets have nearly spherical shapes and are separated by 
a nesting wave- vector Q = (0, 0, 1), which is equivalent to 
(1,0,0) in the Brillouin zone of a body-centered tetrag- 
onal latticed Direct calculation of a static momentum- 
dependent susceptibility 30 also shows a peak at the com- 
mensurate wave-vector Q = (1,0,0). 

A fast decrease of the uniform susceptibility^ 3 , below 
T m as well as suppression of the transition temperature 
T m and the bulk gap Ao by applied magnetic field 3 ^ 3 ^ 3 ^ 
also point to a charge- or a spin-density wave state. For 
the CDW the Zeeman splitting degrades the nesting of 
the Fermi surfaces and reduces a mean-field transition 
temperature 3 ^ in a way, which is analogous to the para- 
magnetic limit effect in superconductors. By contrast, an 
isotropic SDW involves coupling of bands with opposite 
spin and the nesting is not affected by a magnetic field. 
A strong spin-orbit coupling in heavy-fermion materials 
creates momentum dependence of the g- factor. If the 
nesting condition e(k + Q) = — e(k) is satisfied for par- 
ticular sheets of the Fermi surface it is not generally ful- 
filled for the Zeeman shift ^i B g(k+Q)H/2 ^ ^ B g{^)H/2. 
Hence, in metals with strong spin-orbit coupling an SDW 
state is also suppressed by the paramagnetic effect. 

The mean-field theory of a SDW formation in ideally 
nested electron and hole Fermi surfaces 3 ^ resembles to a 
large extent the BCS theory. The relative jump of the 
specific heat at T m — Tsdw is estimated by the BCS 
value AC/C ~ 1.43 , which is compatible with the ex- 
perimental value AC/C « 2.9 once additional strong- 
coupling effects are taken into account. The modulation 
of the spin density at T — are given by 

M Q = ^(4+QT c kT> = HbN A In ^ , (1) 

k 

where Nq is the density of states per one spin direction. 
Estimating Nq ~ n e /ep, we find that ordered moments 
normalized per one electron constitute a small fraction 
of the Bohr magneton ~ A /spi— Such a reduction has 
a transparent physical meaning: only electrons (holes) 
within a thin layer of width 2Aq around the Fermi surface 



participate in the formation of ordered moments. Thus, 
affecting strongly thermodynamic and kinetic properties, 
a weak-coupling SDW order has a small form-factor and 
does not produce significant magnetic Bragg scattering. 
This fact has not been so far appreciated in the literature 
on URu2Si2- 

There are several additional factors, which complicate 
the simple picture drawn above. First, the perfect nesting 
between different bands does not appear in real materials. 
Absence of perfect nesting acts as a depairing effect re- 
ducing gradually both the transition temperature and the 
zero-T gap and enhancing the residual density of states. 
Obviously, such an effect does not change the conclusion 
about a small form-factor, but may reduce the jump in 
the specific heat compared to the BCS value. In order 
to show that partial nesting does not modify the previ- 
ous estimate, we refer to a similar situation in supercon- 
ductors with paramagnetic (depairing) impurities. Using 
the Abrikosov-Gor'kov theory, Skalski et al~ have cal- 
culated the effect of paramagnetic impurities on various 
characteristics of an s-wave superconductor. Their re- 
sults indicate that in a wide range of impurity concentra- 
tion, the jump in the specific heat and the transition tem- 
perature are suppressed at approximately the same rate, 
hence, preserving the BCS estimate for the relative spe- 
cific heat jump. Second, the electron mass enhancement 
in heavy fermion materials (m* /m ~ 25 in URu2Si2) can 
significantly reduce the Fermi energy scale ep. However, 
simultaneously with a mass renormalization, an interac- 
tion with spin fluctuations reduces strongly the spectral 
weight of heavy quasiparticles^S^i adding an extra small 
factor (m/m*) 3 to Eq. Q, which completely compen- 
sates the factor (m*/m) in the density of states and fur- 
ther reduces value of the ordered moments. Finally, ac- 
cording to the band structure calculation ;] 30 ! 33 URu 2 Si2 
is a compensated metal with equal number of electrons 
and holes. The number of carriers in the two bands un- 
dergoing a SDW transition is smaller than one, when 
normalized to the number of U-atoms. This yields an 
extra reduction factor, since the neutron scattering ex- 
periments report the ordered moments normalized per 
one uranium. 

The above arguments can, in our view, convincingly 
explain why small antiferromagnetic Bragg peaks in 
URu2Si2 are consistent with a SDW instability. In the 
following we assume that due to nesting between some 
parts of the Fermi surface in URu2Si2 the commensurate 
SDW state is formed below the critical temperature T m 

and that the SDW amplitude ip ~ Sk( c k+QT c kT) Pl a Y s 
the role of a hidden order parameter in the problem. 



III. PHASE DIAGRAM UNDER PRESSURE 

The microscopic dual modelsiiiS assume that conduc- 
tion electrons are noninteracting and operate, therefore, 
with a single order parameter, which leaves no place for 
the phase diagram with SMAF and LMAF regions. In 



our scenario, temperatures of intrinsic phase transitions 
in itinerant and localized magnetic subsystems are dif- 
ferent functions of P and they may interchange their 
order under pressure. As a result, a line of first-order 
transitions appears naturally between the two ordered 
states, where one order parameter prevails over another, 
see Fig. ^ 

The Landau free-energy for two interacting order pa- 
rameters can be written as 

F = a 1 tp 2 +a 2 m 2 +2-ftpm+Pitp 4 +p 2 m 4 +2p t tp 2 m 2 . (2) 

A special bilinear coupling term is allowed only if the 
two parameters transform according to the same irre- 
ducible representation, otherwise 7 = 0. For nonzero 7, 
the quantities tp and m do not correspond to two differ- 
ent types of symmetry breaking. Rather they describe 
two weakly coupled magnetic subsystems of URu2Si2 in 
a way, which is reminiscent of the Ginzburg-Landau de- 
scription of the multigap superconductivity in MgB2. 31 
The bilinear term corresponds, then, to a polarization of 
local moments by a spin-density wave. 

Generally, in addition to the bilinear term tpm there 
are possible other coupling terms in the Landau func- 
tional: tp 3 m and 1pm 3 . These terms can exist even if tp 
and m transform according to different irreducible repre- 
sentations (though tp has to break the time-reversal sym- 
metry). The tp 3 m term leads, for example, to a small 
antiferromagnetic component in a state with tp ^ 0, even 
if 7 = 0. The induced m component grows in such a case 
asm~ (T m — T) 3 / 2 , while the neutron diffraction exper- 
iments find a standard mean- field exponent 1/2&L This 
observation suggests that the tpm coupling plays a dom- 
inant role and, hence, the phcnomcnological coefficients 
for tp 3 m and tpm 3 terms have the same smallness as 7. 
In such simple linear transformation allows to 

exclude such terms form the Landau functional without 
modifying significantly the physical meaning of tp and m. 

For 7 = the functional has a form commonly 
found in the theory of phase transitions. Assuming that 
only coefficients ax,2 depend on temperature and pres- 
sure, the energy J2J) describes a phase diagram with two 
crossing lines of second-order transitions determined by 
ai t2 (P, T) — 0. The transition line from a paramagnetic 
state T m (P) has a kink at the crossing point. Presence 
and nature of extra transitions in the order state, where 
ai,a 2 < 0, depend on quartic terms. For Pi < \fP1P2 
or for a weak repulsion between tp and m, there are two 
other lines of second-order transitions transition emerg- 
ing from the crossing point. They separate two states 
with pure ordering, i.e., tp ^ 0, m = and tp = 0, 
m 7^ 0, from a mixed phase with tp ^ and m =/= 0. 
Thus, the phase diagram in this case has a tetracritical 
point. For Pi > y 'P\P 2 or for a strong repulsion between 
two components, the mixed phase does not appear. In- 
stead, there is a single line of first-order transitions in the 
P-T plane given by a\P 2 = o^Pi-, which approaches the 
kink (crossing) point from the ordered side, see the left 
panel in Fig. ^ 
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FIG. 1: The phase diagram of the two-order parameter Lan- 
dau functional for 7 = (left panel) and for 7 7^ (right 
panel) . 

In the following we discuss effect of nonzero 7 on the 
two order parameter functional J2J): problem, which, to 
our knowledge, has not been considered so far. Once 
7^0, the two order parameters appear simultaneously 
on a single transition line given by ot\a 2 — 7 2 . The tran- 
sition temperature from a paramagnetic state T m (P) has 
now a smooth pressure dependence and does not exhibit 
a kink. At P = 0, the induced antiferromagnetic compo- 
nent behaves as 

m^-(-f/a 2 )tp (3) 

for I cki j , j'y I <C a 2 - A small coefficient 7/^2 implies 
weak ordered moments, while tp gives rise to a large 
anomaly in the specific heat at T m . We identify this 
phase with a small moment antiferromagnetic (SMAF) 
phase of URu2Si2- 

In order to investigate the possible ordered states and 
phase transitions below T m (P), one has to minimize 
Eq. (0) with respect to both tp and m. This gives a 
system of two coupled cubic equations, which is easily 
solved numerically, but does not allow full analytic so- 
lution. Still, simple analytic arguments can be used to 
prove stability of the first-order transition line for 7^0 
and Pi > y 'PiP 2 . [For Pi < \J P\P 2 , the bilinear term sta- 
bilizes the mixed phase (tp ^ and m ^ 0) right below 
T m (P).} 

Substitution m — > (Pi/P^ 1 ^ 4 m transforms the free- 
energy to a more symmetric form: 

F = a!tp 2 + a 2 m 2 + 2jipm + Pi(tp 2 + m 2 ) 2 + 2p i tp 2 m 2 , 

(4) 

where & 2 = (Pi/P 2 )^ 2 a 2 , 7 = (Pi/p 2 )^ 4 j, and ft = 
{Pi/ P2) 1 / 2 Pi — Pi- In the new notations the condition 
for the absence of the mixed phase at 7 = is Pi > 0, 
while the position of the first-order transitions line in the 
P-T plane is given by ot\ = a 2 . Let us cross from the 
paramagnetic state into the ordered state along this line 
taking a± — a 2 — a. The transformation tp = r\\ — rj 2 , 
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m = r/i + r/2 diagonalizes the quadratic terms in Eq. Q 
yielding 

F = 2(«+7)^+2(«-7) % 2 +4/3 1 (^ +r/ 2 ) 2 +2/ 5 i(7/ 2_ 7? 2 ) 2 _ 

If we assume, for example, 7 < 0, then a second order 
transition takes place at a(T m ) = —7 from a paramag- 
netic state into a state with r/f = —(a + 7)/2(2p\ + 
while 772 = 0. For positive p\, the last term in Eq. © 
disfavors states with rjl 7^ ry|. Therefore, at sufficiently 
low temperature there should be another transition into 
a state with a nonzero 772. The location of such a critical 
point (T c , P c ) is given by 



A A - VP1P2 



(6) 



The ratio of the specific heat jumps at two consecutive 
transitions T m and T c is 



A(C/T) C = A_ 
A(C/T) m 2/3i " 

For a < a(T c ) the two components behave as 



Vl = = ! % = = 

8p\ft 8p\A 



(7) 



(8) 



The relative phase between 771 and 772 is not fixed, though 
solutions ( 1 771 1 , 1 T72 1 ) and (|?7i|, — \i}2\) describe two essen- 
tially different states. Away from the line a\ =0.2, the 
energy (JSJ acquires the extra term 2 (a 1 — 0:2)771 772, which 
immediately lifts the above two-fold degeneracy and se- 
lects either or tt shift between 771 and 772 on the two 
sides of oi = 02- Consequently, the first order transi- 
tion line Tm{P) is stable and its position in the P-T 
plane is given by the same equation as for 7 = 0. How- 
ever, Tm(P) splits from the line of second order phase 
transitions T m (P) and terminates at the critical point 
determined by Eq. ©, see the right panel of Fig.^ 

The two states to the left and to the right from 
Tm(P) are phases with large ipL = \ + 1*72 1 and small 
ms = — 1^2 1 (SMAF) and with small ips — \r)i\ — \r)2\ 
and large = |?7i| + 1 772 1 (LMAF). A relative jump of 
the ordered antiferromagnetic moments across Tm(P) is 
given by 



m L - m s _ |t?2_| 
m L +m s 1 771 1 



o(T c ) 



a + a{T c ) 



(9) 



IV. CRYSTAL-FIELD MODEL FOR INDUCED 
MOMENTS 

The nine-fold degenerate state of U 4+ ions with the to- 
tal angular momentum J = 4 is split by a crystalline elec- 
tric field. Following the previous workf&iLiSi we assume 
that the ground and the first excited levels are singlets 
separated by a crystal field gap A and that the only non- 
vanishing matrix element of the total angular momentum 
is (0|J Z |1) = [i. Working in the basis of the two lowest 
levels, the new pseudo-spin-1/2 operators are defined as 



s*|o) = +i|o), s z 



i) 



(10) 



The nonzero component of the angular momentum oper- 
ator is expressed in terms of pseudo-spin-1/2 operators 
as J z — 2(iS x . Local moments formed by the mixing of 
two crystal field levels have, therefore, a very anisotropic 
nature. They couple only to a z-component of an applied 
field and via Ising-like interaction between different sites. 
Coupling between local moments and itinerant carriers is 
described by an Ising-Kondo interactionii 



C ia a af3 C i[3 ■ 



(11) 



Below the SDW transition, space modulation of the elec- 
tron spin density produces an internal staggered field on 
uranium sites: 



H s (v i )=H s e 1 ^ , H s 



-IM, 



Q 



7/, 



(12) 



The estimate of the staggered field from the experimental 
data on URu2Si2 is given in section V. Transformation 
from Eq. I|llf> to (|12|) corresponds to a mean-field approx- 
imation, which should be sufficient when considering two 
weakly interacting subsystems. 

The total crystal-field Hamiltonian in the presence of 
both staggered and uniform fields applied parallel to the 
crystal z-axis is written in terms of pseudo-spin operators 
as 

H = l,r» i; ' ;S '^ A »< 

(id) « 

-2fij2( H + H ^ Qri ) s ?> (1 3 ) 

i 

where J{ij) is a set of exchange constants between local 
moments on a body centered tetragonal lattice. 



The size of the jump varies continuously along Tm(P) 
and vanishes at P = P c . Note, that the distance between 
the critical line T m (P) and the critical point T c given by 
Eq. © is proportional to 7 and may be quite small. At 
present, the neutron experiments under hydrostatic pres- 
sure has failed to identify the critical end point (T c , P c ) 
on the line of first-order transitions Tm(P)^ We sug- 
gest that specific heat measurements can help to finally 
resolve the phase diagram of URu2Si2- 



A. Zero-field case 

The crystal field model (|13fl in zero applied field has 
been studied by many authors ^^SJU^^ The Hamilto- 
nian 113fl without the last term corresponds to a ubiqui- 
tous Ising model in a transverse field. In order to make 
connection with previous works we briefly list in this sub- 
section the main results on the crystal-field model i|13|) 
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with H = H s = 0. At zero temperature and in the 
large gap limit the system remains in a singlet ground 
state. The excitation spectrum consists of magnetic or 
Van-Vleck excitons, which are bound states of the two 
singlet levels. Energy of Van-Vleck excitons is easily 
found by applying the Holstein-Primakoff representation 
to pseudo-spin-1/2 operators. In the harmonic approxi- 
mation it is suffice to write 



CZ 1 t 



QX 



,t 



The excitation spectrum is given by 



uj k = \/A(A + 2A) 



(14) 



(15) 



where J7k = JV j is a Fourier transform of the 

exchange interaction. The excitation gap is reduced by 
magnetic interactions to 



A g = VA(A - A c ) , A c = 2// 2 |JqI 



(16) 



where the wave-vector Q corresponds to the absolute 
minimum of J^. Let us emphasize here that neither 
A nor A s have any relation to the bulk gap Ao, which 
corresponds to itinerant magnetic subsystem. Quantum 
fluctuations somewhat renormalize the spectrum l|15ll at 
T = and tend to further reduce the critical gap A c 4^£ 
This effect depends, however, on a particular form of J7k, 
and for a three-dimensional system does not exceed 2- 
3%. Below, we neglect such quantum corrections. 

If the crystal field splitting A becomes smaller than A c , 
the system develops a long-range magnetic order with 
a staggered magnetization (J?) = 2fi(Sf) ~ e lCJri . In 
the following we always assume that magnetic ordering 
has a two-sublattice antiferromagnetic structure, that is 
e 2«Qr, = i or 2Q = as in URu 2 Si 2 . In order to describe 
a finite-temperature transition into ordered state one can 
use a simple molecular-field approximation^ 2 ^ For this 
we write 



(Sf) 



(17) 



where the dimensionless staggered magnetization m s is 
determined by a self-consistency equation obtained from 
a single-site solution: 



777 



2^ 2 |JqI 



A 2 + 16/A72 TO 2 



I A 2 + 16 fi 4 J 2 m 2 

:tanh-^ . (18) 

2T v ' 



The transition temperature obtained from the above 
equation is 



A , A c + A 
7^— = In — — . 

TV A c - A 



(19) 



In the molecular-field approximation the transition tem- 
perature vanishes as A — > A c — in agreement with 



Eq. i|16|) . At zero temperature the sublattice magneti- 
zation is 



M. 



sO 



2fj,m s o 




(20) 



whereas near XV the antiferromagnetic moments follow 
the mean-field temperature dependence: 



Mi k, Mi 



T N -T 



T N A, 



- ±A 2 

2 c 



\A 2 



X 



A" 



(21) 



The excitation spectrum in the ordered phase at zero 
temperature is found by introducing a staggered canting 
angle ip for the two sublattices4& In the mean-field ap- 
proximation cos ip — A/A c . After transformation to the 
local (rotating) frame and application of Eq. (|14|) one 
finds: 



w k = 2/j 2 |J7q 




(22) 



The above equation shows that upon approaching the 
Ising limit A -C A c , the dispersion of the longitudinal 
excitations is gradually diminished. For more details and 
discussion see the end of the subsection C. 

Neutron scattering measurements on URu2Si2 yield a 
moderate value of the matrix element of the total angular 
momentum [i ~ 1.2/ib^ A simple explanation of small 
static moments would be, then, to assume that (A c — 
A) <C A. According to Eq. I|19fl such an assumption 
also leads to a small transition temperature compared to 
the crystal-field level splitting TV <C A, which is again 
in agreement with the experimental observation of A ~ 
10meVi& The above straightforward explanation of small 
ordered moments fails, however, to explain a large jump 
of the specific heat at Tm- Indeed, in the molecular-field 
approximation the specific heat jump at the transition 
temperature ifHfl) is 



AC 



2T 



N 



A 2 



dm 2 . 



dT 



Using Eq. (J2TJ we find in the limit A — > A c 

AC _ Ag- A 2 A c + A 
~C ~ 2A 2 n A c - A ' 



(23) 



(24) 



Taking M s q » 0.03/jb, which implies that (A c — A)/A c sw 
3 x 10~ 4 , we find for the specific heat jump AC/C w 
2.7 x 1CV 3 . Such a jump is three orders of magnitude 
smaller than the experimentally measured jump at the 
17.5 K transition^ Corrections to the molecular-field 
approximation^^ do not significantly modify the jump 
AC. Consequently, it has been concluded that sponta- 
neous ordering of local moments on uranium sites can- 
not explain the phenomenology of the antiferromagnetic 
transition in URU2S12. In the next sections we shall con- 
sider the model (|13|) in the regime of induced local mo- 
ments, that is A > A c = 2/j 2 |Jq| and H s ^ 0. 
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B. Finite Fields: Mean-field approximation 

The mean-field ansatz for a sublattice magnetization in 
the presence of both uniform H and staggered H s mag- 
netic fields is given by 



(Sf) = m s e 1 ^ + mo 



(25) 



For a single spin, the mean-field Hamiltonian takes the 
following form 

H M F = -AS? -S*[(2nH s +4n 2 \J Q \m s y Qr > 

+ 2 f iH - 4fi 2 J mo] , (26) 

where j7o = J7k=o- Calculating an equilibrium magneti- 
zation we find for two sublattices: 



m s ± mo 



D± 



: tanh 



A 2 + 4D£ 2T 

2/ 



D ± = ,i(ff s ±H) + 2^(\J Q \m s T J m ) . (27) 

Below, we focus on the case A > A c = 2/x 2 |J7q|, when 
there is no magnetic ordering in the subsystem of local 
moments down to T — in the absence of both external 
and internal fields. For weak staggered field, linearization 
of Eq. (|2*7|) in H s and m s at H = yields 



/ii? s tanh(A/2T) 
A - A c tanh(A/2T) 



(28) 



In accordance with the phenomcnological formula of 
section II, weak local moments are proportional to the 
primary (SDW) order parameter. At zero temperature 
the dimensional staggered moments are 



M s0 = 2/x- 



fiH s 



(29) 



The above equation allows to estimate the staggered field 
in URu2Si2 from the available experimental data, see sec- 
tion V. 

The effect of a uniform field on induced local moments 
is considered, for simplicity, for T — only. In this case 
expansion of Eq. i|27|) to linear order in m s and H s yields 



A(l + Vff 2 /A 2 ) 3 / 2 - A c 



(30) 



where an effective field H = H — inJ^ra® is determined 
self-consistently from 



H = H 



2^ 2 J H 



(A 2 + ViJ 2 )V2 



(31) 



According to arguments of section III, the SDW con- 
tribution to the magnetic Bragg peaks is negligible due 
to a small form-factor. Then, the measured intensity 
of Bragg reflections is proportional to to 2 . Suppression 



of an SDW order parameter with an external field can 
be described by a simple formula ip 2 cx (1 — H 2 /H 2 ) or 
H 2 = H%(1- H 2 /H 2 ), where H c ~ 40 T is a meta- 
magnetic field in URu2Si2^ The magnetic Bragg peak 
intensity is 



y 2 H 2 (l-H 2 /H 2 ) 
[A(l + 4 M 2 iJ 2 /A 2 ) 3 / 2 - A c ] s 



(32) 



This zero-temperature result should be compared to the 
analogous formulas valid near T m , which have been de- 
rived in the previous works 20-24 from the Landau free- 
energy functional. Though different in the details, the 
two limits exhibit an inflection point in Iq(H). The stag- 
gered magnetization remains finite until H c , when the 
primary order parameter is suppressed to zero. The or- 
dered moments are, however, substantially reduced com- 
pared to its zero-field value at significantly smaller mag- 
netic field. Indeed, expanding Eqs. IpTT)) and (|3^|l to the 
first order in H 2 we obtain 



m 2 a (H) 

TO 2 (0) 



1 

m 



12/i 2 A 



(A-A c )(A+2/i 2 Jo) 2 



(33) 



For the completeness we also note that the ferromag- 
netic component of the induced magnetic moments is 
given by 



toq 



A + 2 M 2 Jo 



(34) 



for fields smaller than H* — (A + 2/^ 2 j7o)/2/i. Above this 
field the ferromagnetic component remains constant un- 
til a metamagnetic transition related to a crossing with 
higher energy crystal-field levels. The uniform compo- 
nent of the induced moments too should be measurable 
from a magnetic contribution to the nuclear Bragg peaks. 



C. Finite Fields: Energy Spectrum 

We start with the case H ^ 0, H s = 0, since an ef- 
fective H s in URu2Si2 should be quite small. Partial po- 
larization of magnetic moments (pseudo-spins) along z 
(x) axis is taken into account by rotation of pseudo-spins 
from a laboratory frame to a local (primed) frame: 



Sf = Sf cos p + sin p , 
Sf = -Sf' sin p> + Sf' cos ip . 



(35) 



In the transformed frame and omitting primes the Hamil- 
tonian (|13fl takes the following form: 



H = 4^i 2 V J{ij) Sf Sf cos 2 '■y + SfSUh?y 



(hi) 



+(SfS] +SfSf)sincp cos cp 
- ^2 [( A cos f + 2 M# sin S i 

i 

+ (2^H cos tp - A sin p) Sf 



(36) 
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The boson representation (|14fl of the pseudo-spin opera- 
tors is applied to the above Hamiltonian and the rotation 
angle ip is determined from the condition of vanishing lin- 
ear terms in a,- and aj: 



A tan ip + 2/i To sin ip = 2/iH 



(37) 



At small fields ip w 2^H/(A + 2/i 2 J~o)- 

The harmonic part of the Hamiltonian l|3tj|> after the 
Fourier transformation becomes 



H 2 = a^aktA cos ip + 2fiH sin ip — 2/i 2 sin 2 y>j7o 
1 



k 



/i 2 cosVjk] + ^ // 2 cos 2 ( ) c 1 7k(aka-k + akaL k )- (38) 

The k-independent term is simplified with the help of 
Eq. (j37(l to A/ cos<^ and after applying the Bogoliubov 
transformation we obtain the following field dependence 
of the exciton spectrum: 
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— — + 2fj, 2 J k A cos ip 

COS z (p 



(39) 



The gap at k = Q increases quadratically with magnetic 
field: 



i 2 / TT\ ~ , A / A A^ 



(A + 2n 2 Jo) : 



■A(2A + A C ) . (40) 



The parabolic law for the field dependence of the 
gap has recently been measured in neutron scattering 
experiments^ For an arbitrary wave- vector the field de- 
pendence of the exciton energy is 



ujI(H) « c 2 (0) 



4/i 2 ff 2 A 

(A + 2^j y 



[A - f i 2 J k ] . (41) 



The field dependence changes its sign, i.e., the energy 
starts to decrease with magnetic field, for the wave- 
vectors such that /i 2 7k > A. In terms of zero field fre- 
quencies this is equivalent to Wk > \/3A. In the region in 
the Brillouin zone where Uk ~ \/3A the field dependence 
of the spectrum becomes vanishingly small. Experimen- 
tally, a drastic change in the field dependence has been 
observed between k = Q = (1,0,0) and k = (1.4,0,0). 
The present theory explains a qualitative difference in 
the field response of the two types of excitons. Further 
inelastic neutron measurements on URu2Si2 should allow 
a detailed comparison with our theory and extraction of 
microscopic parameters from experimental data. 

If both staggered and uniform fields are present, the 
derivation of the spectrum becomes a bit more com- 
plicated. One has to explicitly introduce two types of 
bosons <2j and bi for two antifcrromagnetic sublattices 
and to calculate spectrum in the magnetic Brillouin zone, 
which is twice smaller than an original lattice Brillouin 
zone used above. The transformation to local (primed) 
axes from the laboratory frame is given by 

Sf = Sf cos (p, + S z 'e lClr ' simp, , 

S* = -Sf' e iQr * sin (pt + S? cos w , (42) 



where two angles ipi = <pi, for e lC * Ti = 1 and ipi = (p2, 
for e lC * ri — — 1 describe different response of the two 
sublattices. The angles are determined by 

Atan^i + 2fi 2 Ji sintpi = 2fi(H s + H) + 2[i 2 J 2 smtp 2l 
A tan tp2 + 2p?J\ sin</? 2 = 2fj,(H s —H) + 2^ 2 J 2 sin^i. 

here we defined separate summation of exchange con- 
stants over the same J\ — J^i jeA 3^3) anc ^ the different 
sublattices J 2 = J2ieA,jeB ^fe')- 

Performing the same steps as in the case of H s = we 
in the end find 



.±2 



1 



I 1 i 2 \ 

;Kk + w 2k ) 



Ml \ ( W lk - W 2k) 2 + 4 M 4 ^2k A2 COS ¥>1 C0S 



A 2 



Mk cos 2 tp 12 



2fi 2 Ji k A cos ip lj2 ■ (43) 



Here, the Fourier transforms are given by J7ik = 

£ ijs ':'./i'' kr and ^k = V, B^foV" ■ 
The characteristic feature of this spectrum is a small 
jump between two branches of excitations u>£ and ui^ 
at the magnetic Brillouin zone boundary, where 72k = 0. 
In URu 2 Si2 {H s ^ 0) such a jump is induced by external 
magnetic field H ~ H s and becomes negligible again for 
H 3> H s , when the above expression Eq. is valid. 

Finally, let us comment on the longitudinal polariza- 
tion of magnetic excitons detected experimentally 8 In 
the harmonic approximation, the dynamic structure fac- 
tor S zz (r, t) = (J?(t)J? +r (t+T)) is expressed via pseudo- 
spin operators as 



S"(t, r) « 4 M 2 cos 2 <p(Sf(t)Sf +r (t + r)) 



(44) 



In a weakly polarized Van-Vleck paramagnet for A ~ A c 
one has costpi.2 ~ 1. Therefore, transverse 'spin-wave- 
type modes' in the pseudo-spin representation corre- 
spond to longitudinal polarization of magnetic excitons. 
The omitted terms in Eq. (Hjl. such as (S*S* +T ) and 
(SfS z +r ), describe a higher energy two-magnon contin- 
uum and its interaction with single-particle states. These 
terms do not modify the conclusion about longitudinal 
polarization of single-particle excitations. In the oppo- 
site limit A < A c in a phase with large antiferromag- 
netic moments cos(p±. 2 — > 0, and the longitudinal dy- 
namic structure factor S zz (q, uj) does not have contri- 
bution from magnetic excitons. This explains why the 
neutron scattering measurements 2 ^ failed to observe the 
magnetic excitations above the first-order transition at 
Pm = 5 kbar. 



V. COMPARISON WITH EXPERIMENT 

Theoretical predictions of the above section can be di- 
rectly compared with the available experimental data. 
Specifically, let us consider the field dependence of the 
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Field [T] 



FIG. 2: Field dependence of the intensity of the antiferromag- 
netic Bragg peak at Q = (1, 0, 0). Points are the experimental 
datai— Lines are theoretical curves described by Eq. llo21i with 
A = 6 meV (full line) and A = 3 meV (dashed line). 



intensity of the magnetic Bragg peak. The metamagnetic 
transition in URu2Si2 at H c ~ 40 T can be chosen as the 
critical field for the spin-density wave. The field depen- 
dence of the Bragg peak intensity described by Eqs. i|3T|) 
and loL'l) is, then, determined by three microscopic pa- 
rameters: A, A c = 2/i 2 | 1 /q|, and 2/i 2 j7o- The last two 
parameters are fixed by using the experimental data^a for 
the excitation gap A g = ^/A(A — A c ) w 1.59 meV and 
its dependence on applied magnetic field. In this way we 
are left with only one free parameter: the crystal-field 
splitting A. 

The two theoretical curves for A = 3 meV (A c = 
2.2 meV and 2fi 2 J = 1.3meV) and A = 6 meV 
(A c = 5.57 meV and 2^x 2 J = 2.9meV) are presented 
in Fig. together with the experimental datai^i Both 
curves exhibit behavior with an inflection point. The 
larger value of the gap gives better agreement with the 
experimental results. For A = 6 meV the top of the exci- 
ton band at H = calculated from Eq. I|15l) corresponds 
to ljq « 7.3 meV. The internal staggered field estimated 
from Eq. I|29|l in this case is H s as 0.08 T, which is indeed 
much smaller than applied magnetic fields and justifies 
the used approximations. 

If we further increase A, the theoretical dependence 
for Iq(H) with the above two constraints practically 
saturates at the position given by A = 6 meV curve. 
Thus, while we can definitely exclude smaller values 
A < 6 meV for the crystal-field level splitting, the larger 
values A > 6 meV are equally possible. For example, 
for A — 10 meV, which has been suggested on the ba- 
sis of the early neutron scattering measurements^ the 
parameters obtained from the fits are A c = 9.7 meV, 
2/i 2 Jo = 4.9meV, and ujo « 12.2 meV. At present, there 
is no agreement on the value of uiq between the two 
groups of inelastic neutron measurements*^^ 8 . Addi- 
tional more precise neutron scattering investigation of 



URri2Si2 should greatly help to settle this dispute. 



VI. CONCLUSIONS 

We have presented a theoretical model to describe un- 
usual magnetism in URu2Si2 below T m — 17.5 K, which 
combines tiny ordered moments /i ~ 0.03/^s with a large 
specific heat anomaly at the transition point. At ambient 
pressure, the transition is driven by an SDW instability 
in the itinerant subsystem, which also induces weak or- 
dering of local moments on uranium sites. We argue that 
such a low-T c spin-density wave has a small form-factor 
and does not contribute significantly to the neutron scat- 
tering, which essentially probe the localized subsystem. 
Phenomenologically, the phase diagram of URu2Si2 is de- 
scribed by the two order parameter functional (J2J , which 
is consistent with the first-order transition into a state 
with large antiferromagnetic moments. The microscopic 
origin of a strong repulsion between two order parameters 
of the same symmetry needs further clarification. In our 
view such a behavior may result from a strong renormal- 
ization of the RKKY type interaction between the local 
moments by a rather large SDW gap, which opens over 
a half of the Fermi surface. 

In our discussion we have assumed following the previ- 
ous works&iii^ that uranium ions are in the 3 if 4 ground 
state. The validity of such an assumption needs fur- 
ther clarification. Also, the high field behavior with an 
itinerant electron metamagnetic transition 47 would be 
an interesting test for the dual model and its extension 
suggested in in the present work. Another open ques- 
tion is temperature evolution of the crystal-field excita- 
tions. We believe that the experimentally observed dis- 
appearance of the magnetic excitons above the transition 
temperatur e) 8 ! 24 is largely related to the closeness of two 
energy scales: T m = 17.5 K and A g = 1.6 meV. 

The analysis presented in section IVa may be also rel- 
evant to UPta, another heavy-fermion compound with 
tiny antiferromagnetic moments, for review see Ref. l49l 
This material does not have apparent anomalies in ther- 
modynamic and kinetic properties at T m «5K, though 
the neutron diffraction experiments have detected small 
antiferromagnetic moments \x ~ 0.02/1^. In a possible 
scenario for UPt3, there is no a SDW instability in the 
conduction subsystem. The phase transition is driven 
by the RKKY or the superexchange interaction between 
local moments, which only slightly exceeds the criti- 
cal value for zero-temperature antiferromagnetic ordering 
determined by a crystal-field level splitting. While the 
crystal level structure is not precisely known for UPt3, 
the estimates for the specific heat anomaly given in the 
end of section IVa should be generally valid. Thus, the 
small ordered antiferromagnetic moments can be recon- 
ciled with the absence of large anomalies at the transition 
point. The pressure effect on antiferromagnetic order- 
ing in UPt3 also agrees with Eqs. (|19f) and l|20|) . which 
predict a much faster square-root suppression of zero- 
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temperature moments compared to a slow logarithmic 
decrease of the transition temperature. 
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